#########################################################################################################
######################################## STATA REPLICATION CODE #########################################
################################# ANGELO VITO PANARO & ANDREA VACCARO ###################################
### INCOME INEQUALITY IN AUTHORITARIAN REGIMES: THE ROLE OF POLITICAL INSTITUTIONS AND STATE CAPACITY ###
#########################################################################################################
#########################################################################################################

## Load dataset
use "Dataset"

## Set data as panel data
xtset GWn year

## Drop if country has not been classified as autocratic for more than 10 years from 1972 to 2014
keep if number_autocr > 10

## Figure 1
graph box Capacity if ht_regtype2 != 5, over(ht_regtype2)


------------

## Load dataset
use "Dataset"

## Set data as panel data
xtset GWn year

## Figure 2
egen avg_ineq_mon = mean(gini_disp) if ht_regtype2 == 1, by(year)
egen avg_ineq_mil = mean(gini_disp) if ht_regtype2 == 2, by(year)
egen avg_ineq_onep = mean(gini_disp) if ht_regtype2 == 3, by(year)
egen avg_ineq_multip = mean(gini_disp) if ht_regtype2 == 4, by(year)
twoway line avg_ineq_mon year if year <= 2014, sort ylabel(, nogrid) name(ineq_mon1)
twoway line avg_ineq_mil year if year <= 2014, sort ylabel(, nogrid) name(ineq_mil1)
twoway line avg_ineq_onep year if year <= 2014, sort ylabel(, nogrid) name(ineq_onep1)
twoway line avg_ineq_multip year if year <= 2014, sort ylabel(, nogrid) name(ineq_multip1)
gr combine ineq_mon1 ineq_mil1 ineq_onep1 ineq_multip1, rows(2) cols(2) ycommon
gr drop ineq_mon1 ineq_mil1 ineq_onep1 ineq_multip1


------------

## Load dataset
use "Dataset"

## Set data as panel data
xtset GWn year

## Figure 3
keep if ht_democr == 0
collapse (firstnm) vdem_country_name (mean) gini_disp Capacity, by(GWn)
pwcorr gini_disp Capacity
twoway (scatter gini_disp Capacity, mlabel(vdem_country_name)) lfit  gini_disp Capacity


------------

## Load dataset
use "Dataset"

## Set data as panel data
xtset GWn year

## Drop if country has not been classified as autocratic for more than 10 years from 1972 to 2014
keep if number_autocr > 10

## Figure 4
set level 90
xtpcse gini_disp l.b4.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m4
coefplot (m4), keep(1L.ht_regtype2 2L.ht_regtype2 3L.ht_regtype2) yline(0) msymbol(o) mcolor(gs8) msize(vhuge) ciopts(color(gs8)) vert coeflabel(1L.ht_regtype2 = "Monarchy" 2L.ht_regtype2 = "Military" 3L.ht_regtype2 = "One-party") ylabel(-.1(.05).1) legend(off) name(h1)
xtpcse gini_disp l.b3.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m4
coefplot (m4), keep(1L.ht_regtype2 2L.ht_regtype2 4L.ht_regtype2) yline(0) msymbol(o) mcolor(gs8) msize(vhuge) ciopts(color(gs8)) vert coeflabel(1L.ht_regtype2 = "Monarchy" 2L.ht_regtype2 = "Military" 4L.ht_regtype2 = "Multiparty") ylabel(-.1(.05).1) legend(off) name(h2)
xtpcse gini_disp l.b2.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m4
coefplot (m4), keep(1L.ht_regtype2 3L.ht_regtype2 4L.ht_regtype2) yline(0) msymbol(o) mcolor(gs8) msize(vhuge) ciopts(color(gs8)) vert coeflabel(1L.ht_regtype2 = "Monarchy" 3L.ht_regtype2 = "One-party" 4L.ht_regtype2 = "Multiparty") ylabel(-.1(.05).1) legend(off) name(h3)
xtpcse gini_disp l.b1.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m4
coefplot (m4), keep(2L.ht_regtype2 3L.ht_regtype2 4L.ht_regtype2) yline(0) msymbol(o) mcolor(gs8) msize(vhuge) ciopts(color(gs8)) vert coeflabel(2L.ht_regtype2 = "Military" 3L.ht_regtype2 = "One-party" 4L.ht_regtype2 = "Multiparty") ylabel(-.1(.05).1) legend(off) name(h4)
gr combine h1 h2 h3 h4, rows(2) cols(2) ycommon l1title(Slope coefficients)
gr drop h1 h2 h3 h4


------------

## Table 2
# Multiparty
xtpcse gini_disp l.b4.ht_regtype2 l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l(1/2).gini_disp, p
estimates store m1
xtpcse gini_disp l.b4.ht_regtype2 l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m2
xtpcse gini_disp l.b4.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m3
estout m1 m2 m3, cells(b(star fmt(3)) se(par fmt(3))) legend label varlabels(_cons Constant) stats(r2 N N_g, fmt(3 0)) starlevels( * 0.10 ** 0.05 *** 0.010) drop(*GWn *decadevar)
# One-party
xtpcse gini_disp l.b3.ht_regtype2 l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l(1/2).gini_disp, p
estimates store m1
xtpcse gini_disp l.b3.ht_regtype2 l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m2
xtpcse gini_disp l.b3.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m3
estout m1 m2 m3, cells(b(star fmt(3)) se(par fmt(3))) legend label varlabels(_cons Constant) stats(r2 N N_g, fmt(3 0)) starlevels( * 0.10 ** 0.05 *** 0.010) drop(*GWn *decadevar)
# Military
xtpcse gini_disp l.b2.ht_regtype2 l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l(1/2).gini_disp, p
estimates store m1
xtpcse gini_disp l.b2.ht_regtype2 l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m2
xtpcse gini_disp l.b2.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m3
estout m1 m2 m3, cells(b(star fmt(3)) se(par fmt(3))) legend label varlabels(_cons Constant) stats(r2 N N_g, fmt(3 0)) starlevels( * 0.10 ** 0.05 *** 0.010) drop(*GWn *decadevar)
# Monarchy
xtpcse gini_disp l.b1.ht_regtype2 l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l(1/2).gini_disp, p
estimates store m1
xtpcse gini_disp l.b1.ht_regtype2 l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m2
xtpcse gini_disp l.b1.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m3
estout m1 m2 m3, cells(b(star fmt(3)) se(par fmt(3))) legend label varlabels(_cons Constant) stats(r2 N N_g, fmt(3 0)) starlevels( * 0.10 ** 0.05 *** 0.010) drop(*GWn *decadevar)


------------

##### SUPPLEMENTARY MATERIAL #####

## Table S1
xtpcse gini_disp l.b4.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1).gini_disp i.GWn, p
estimates store m1
xtpcse gini_disp l.b4.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m2
xtpcse gini_disp l.b4.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/3).gini_disp i.GWn, p
estimates store m3
xtpcse gini_disp l.b4.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/4).gini_disp i.GWn, p
estimates store m4
xtpcse gini_disp l.b4.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/5).gini_disp i.GWn, p
estimates store m5
estout m1 m2 m3 m4 m5, cells(b(star fmt(3)) se(par fmt(3))) legend label varlabels(_cons Constant) stats(r2 N N_g, fmt(3 0)) starlevels( * 0.10 ** 0.05 *** 0.010) keep(*gini_disp)


------------

## Table S2
xtpcse gini_disp l.b4.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1).gini_disp i.GWn, p
estimates store m1
xtpcse gini_disp l.b4.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m2
xtpcse gini_disp l.b4.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/3).gini_disp i.GWn, p
estimates store m3
estout m1 m2 m3, cells(b(star fmt(3)) se(par fmt(3))) legend label varlabels(_cons Constant) stats(r2 N N_g, fmt(3 0)) starlevels( * 0.10 ** 0.05 *** 0.010) drop(*GWn *decadevar)
xtpcse gini_disp l.b3.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1).gini_disp i.GWn, p
estimates store m1
xtpcse gini_disp l.b3.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m2
xtpcse gini_disp l.b3.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/3).gini_disp i.GWn, p
estimates store m3
estout m1 m2 m3, cells(b(star fmt(3)) se(par fmt(3))) legend label varlabels(_cons Constant) stats(r2 N N_g, fmt(3 0)) starlevels( * 0.10 ** 0.05 *** 0.010) drop(*GWn *decadevar)
xtpcse gini_disp l.b2.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1).gini_disp i.GWn, p
estimates store m1
xtpcse gini_disp l.b2.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m2
xtpcse gini_disp l.b2.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/3).gini_disp i.GWn, p
estimates store m3
estout m1 m2 m3, cells(b(star fmt(3)) se(par fmt(3))) legend label varlabels(_cons Constant) stats(r2 N N_g, fmt(3 0)) starlevels( * 0.10 ** 0.05 *** 0.010) drop(*GWn *decadevar)
xtpcse gini_disp l.b1.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1).gini_disp i.GWn, p
estimates store m1
xtpcse gini_disp l.b1.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m2
xtpcse gini_disp l.b1.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/3).gini_disp i.GWn, p
estimates store m3
estout m1 m2 m3, cells(b(star fmt(3)) se(par fmt(3))) legend label varlabels(_cons Constant) stats(r2 N N_g, fmt(3 0)) starlevels( * 0.10 ** 0.05 *** 0.010) drop(*GWn *decadevar)


------------

# Load dataset
use "Dataset"

# Set data as panel data
xtset GWn year

## Table S3
xtpcse gini_disp l.b4.ht_regtype2 l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l(1/2).gini_disp, p
estimates store m1
xtpcse gini_disp l.b4.ht_regtype2 l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m2
xtpcse gini_disp l.b4.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m3
estout m1 m2 m3, cells(b(star fmt(3)) se(par fmt(3))) legend label varlabels(_cons Constant) stats(r2 N N_g, fmt(3 0)) starlevels( * 0.10 ** 0.05 *** 0.010) drop(*GWn *decadevar)
xtpcse gini_disp l.b3.ht_regtype2 l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l(1/2).gini_disp, p
estimates store m1
xtpcse gini_disp l.b3.ht_regtype2 l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m2
xtpcse gini_disp l.b3.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m3
estout m1 m2 m3, cells(b(star fmt(3)) se(par fmt(3))) legend label varlabels(_cons Constant) stats(r2 N N_g, fmt(3 0)) starlevels( * 0.10 ** 0.05 *** 0.010) drop(*GWn *decadevar)
xtpcse gini_disp l.b2.ht_regtype2 l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l(1/2).gini_disp, p
estimates store m1
xtpcse gini_disp l.b2.ht_regtype2 l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m2
xtpcse gini_disp l.b2.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m3
estout m1 m2 m3, cells(b(star fmt(3)) se(par fmt(3))) legend label varlabels(_cons Constant) stats(r2 N N_g, fmt(3 0)) starlevels( * 0.10 ** 0.05 *** 0.010) drop(*GWn *decadevar)
xtpcse gini_disp l.b1.ht_regtype2 l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l(1/2).gini_disp, p
estimates store m1
xtpcse gini_disp l.b1.ht_regtype2 l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m2
xtpcse gini_disp l.b1.ht_regtype2 l.Capacity l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m3
estout m1 m2 m3, cells(b(star fmt(3)) se(par fmt(3))) legend label varlabels(_cons Constant) stats(r2 N N_g, fmt(3 0)) starlevels( * 0.10 ** 0.05 *** 0.010) drop(*GWn *decadevar)


------------

# Drop if country has not been classified as autocratic for more than 10 years from 1972 to 2014
keep if number_autocr > 10

## Table S4
xtpcse gini_disp l.b4.ht_regtype2 l.v2xcl_prpty l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m1
estout m1, cells(b(star fmt(3)) se(par fmt(3))) legend label varlabels(_cons Constant) stats(r2 N N_g, fmt(3 0)) starlevels( * 0.10 ** 0.05 *** 0.010) drop(*GWn *decadevar)
xtpcse gini_disp l.b3.ht_regtype2 l.v2xcl_prpty l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m1
estout m1, cells(b(star fmt(3)) se(par fmt(3))) legend label varlabels(_cons Constant) stats(r2 N N_g, fmt(3 0)) starlevels( * 0.10 ** 0.05 *** 0.010) drop(*GWn *decadevar)
xtpcse gini_disp l.b2.ht_regtype2 l.v2xcl_prpty l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m1
estout m1, cells(b(star fmt(3)) se(par fmt(3))) legend label varlabels(_cons Constant) stats(r2 N N_g, fmt(3 0)) starlevels( * 0.10 ** 0.05 *** 0.010) drop(*GWn *decadevar)
xtpcse gini_disp l.b1.ht_regtype2 l.v2xcl_prpty l.ln_wdi_gdppc2010 l.sq_ln_wdi_gdppc2010 l.ln_wdi_totalpop l.wdi_trade l.wdi_natresrents l.v2pepwrses l.b2.dpi_execrlc i.b1.decadevar l(1/2).gini_disp i.GWn, p
estimates store m1
estout m1, cells(b(star fmt(3)) se(par fmt(3))) legend label varlabels(_cons Constant) stats(r2 N N_g, fmt(3 0)) starlevels( * 0.10 ** 0.05 *** 0.010) drop(*GWn *decadevar)


